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ABSTRACT 

Rogue planets have been ejected from their planetary system. We investigate the possibility that a 
rogue planet could maintain a liquid ocean under layers of thermally-insulating water ice and frozen 
gas as a result of geothermal heat flux. We find that a rogue planet of Earth-like composition and age 
could maintain a subglacial liquid ocean if it were faS.S times more massive than Earth, corresponding 
to ~8 km of ice. Suppression of the melting point by contaminants, a layer of frozen gas, or a larger 
complement of water could significantly reduce the planetary mass that is required to maintain a 
liquid ocean. Such a planet could be detected from reflected solar radiation, and its thermal emission 
could be characterized in the far-IR if it were to pass within 0(1000) AU of Earth. 

Subject headings: astrobiology - conduction - convection - planetary systems - planets and satellites: 
surfaces 
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1. INTRODUCTION 

As a planetary system forms, some planets or plan- 
etesimals, referred to as "rogue" planets, can enter hy- 
perbolic orbits and be ejected from the system as a re- 
sult of gravitationa l interactions with gas giant plan- 
ets (jLissauerl Il987f l. Furthermore, interaction with 
passing stars can eject planets from mature systems 
(jLaughfln fc Adamslf2000t) . The ability of a rogue planet 
to support life is of interest as a sort of pathological ex- 
ample of planetary habitability, because such a planet 
could potentiall y represent a viable optio n for interstel- 
lar panspermia (jPurand-Manterolal 120101 ) . and because 
such a planet could be the closest source of extrasolar 
life for exploration by humanity in the distant future. 
Since some sort of starting point is required to discuss 
the issue, a planet is often deflned as habitable if it ca n 
sustain liquid wa ter at its surface (|Kasting et al.l [19931 ). 
iStevensonI ()1999[ ) argued that if a rogue planet had an 
extremely high-pressure hydrogen atmosphere, pressure- 
broadening of far-infrared absorption by molecular hy- 
drogen could support liquid water on the planet's surface 
as a result of the geothermal heat flux alone, making the 
planet potentially habitable. iDebes fc SigurdssonI (|2007f ) 
showed that terrestrial planets can be ejected with moons 
and that the resulting tidal dissipation could increase the 
geothermal heat flux by up to two orders of magnitude 
for 0(10* yr). 

Subglacial liquid water oceans sustained by internal 
heat flux on icy bodies represent an alternative type 
of habitat. It is well known that subglacial oceans are 
possible on moons around giant plane ts and on trans- 
Nept unian objects in the solar system (jHussmann et al.l 
[20061). as weU as water- rich exoplanets in distant orbits 
(jEhrenreich et a l. 200d: iFu et"al]|2010[) . A possible ter- 
restrial analog is Lake Vostok, a «125 m deep lake which 
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is sustained by geothermal heat flu x under fti4 km of ice 
on An tarctica (jKapitsa et al.l [19961 ). iLaughlin fc Adain^ 
([20M l have even argued that a terrestrial rogue planet, 
not attached to any star and receiving negligible energy 
at its surface, could sustain a subglacial liquid ocean if 
it had a thick enough ice layer. We wish to consider this 
point in more depth, including issues such as the poten- 
tial for solid-state convection of ice, the potential effect 
of a thermally insulating frozen gas layer from outgassing 
of the mantle on an Earth-like rogue planet, the effects 
of melting point suppression due to contaminants, and 
observational prospects. By Earth-like, we mean specif- 
ically within an order of magnitude in mass and water 
complement, similar in composition of radionuclides in 
the mantle, and of similar age. A subglacial ocean on 
a rogue planet is interesting because it could serve as 
a habitat for life which could, for example, survive by 
exploiting chemical energy of rock that is continually ex- 
posed by an active mantle. We will refer to a rogue planet 
harboring a subglacial ocean as a Steppenwolf planet, 
since any life in this strange habitat would exist like a 
lone wolf wandering the galactic steppe. 

We can imagine that the ice layer on top of an ocean 
on a Steppenwolf planet will grow until either it reaches 
a steady state with the ice bottom at the melting point, 
or all available water freezes. Geothermal heat from the 
interior of the Steppenwolf planet will be carried through 
the ice layer by conduction, and potentially by convec- 
tion in the lower, warmer, and less viscous portion of the 
ice layer. Since convection transports heat much more ef- 
flciently than conduction, the steady-state ice thickness 
will be much larger if convection occurs, making it harder 
to maintain a subglacial ocean. 

Here we will calculate steady-state ice thicknesses when 
there is conduction only and when there is convection in 
the lower portion of the ice, and make the conservative 
assumption that the thicker solution is valid. We must 
acknowledge, however, that it is very difficult to establish 
definitively whether convection would occur, and the re- 
sulting ice thickness if it were to occur, without detailed 
knowledge of conditions in and rn icroscale composition 
of the ice (jBarr fc Showmaiill2009[ ). More generally, we 
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will make many simplifications, including considering the 
question within the framework of a one-dimensional (ver- 
tical) model, since our primary objective is to establish 
whether or not a Steppcnwolf planet is feasible. 

2. GEOPHYSICAL CONSIDERATIONS 

First we calculate the conductive steady-state thick- 
ness, -ffcond- Above «10 K, the temperature de- 
pendence of the thermal conductivity of water ice 
is well- approximated by k{T) = AT~^, where T 
is the temperature in Kelv in and A=651 W m~^ 
(jPetrenko fc WhitworthI |2002| ). Dimensional analysis 
shows that thermal steady state is reached in ^ 10^ 
years, much shorter than the timescale of decay of the 
geothermal heat flux. Geothermal heat flux through the 
shell will be constant at steady state, since no heat is 
produced within the ice, as would occur by tidal heating 
of a frozen moon. Since the Steppenwolf planets we con- 
sider would be much larger and drier than the icy moons 
on which subglacial oceans are typically studied, we can 
assume that the ice thickness is much less than the plan- 
etary radius, yielding an exponential temperature profile 
through the ice and steady-state thickness, so that 



i^cond = ^ log 1^ ^ 



(1) 



where Th is the temperature at the ice-water interface 
(the melting temperature), Tq is the temperature at the 
top of the ice, and F is the geothermal heat flux. 

Decay of radioactive elements in Earth's interior and 
primordial heat remaining from Earth's formation lead to 
an average geothermal heat fl ux emanating from E arth's 
surface of Fe =0.087 W m-2 (jPollack et al.|[T99l . This 
heat flux decays with time such that Earth's geother- 
mal heat flux may have been roug hly twice its present 
value 3 Gyr ago (jTurcottd Il980f ). In order to con- 
sider Steppenwolf planets of different sizes, we use the 
radius-m ass scaling R oc M "" for super-Earths with 
V = 0.27 (jValencia et al.|[2006[ ). Heuristically, this yields 
a geothermal heat flux that scales as (M/Af^)^"^", or 
roughly as the square root of the mass. 

The pressure at the bottom of the ice layer is « 
9 MPa for each kilometer of ice, scaling with mass as 
(M/M©)i-2^. The mehing point of pure ice is 250-270 K 
at pr essures less than 620 MPa (jChoukroun fc Grassetl 
|2007| ) , although contaminants to pure ice such as chloride 
salts or ammon ia could suppress the melting point by 
» 50 - -100 K (lKargellll99U [19921: iF brtes fc Choukrounl 
120101: iGrasset fc Sotinlll996( ). The steady-state thickness 
in the conductive regime is only logarithmically sensi- 
tive to Th, so we will take Th = 260 K in the estimates 
here, except that we will consider the eutectic point of 
an ammonia- water ice mixture at 176 K to demonstrate 
the impact of contaminants. In steady-state, the tem- 
perature at the surface of a Steppenwolf planet (Tg), i.e., 
the top of the ice or frozen gas layer, will be set by a bal- 
ance between thermal emission and geothermal heat flux, 
F = aT^, where a is the Stefan-B oltzmann constant . 
Astrophysical ra diation backgrounds (jMathis et al.lll983l : 
IDole et al.l[2006h are negligible. 

Any gas present in the atmosphere at planetary ejec- 
tion or outgassed subsequently by geological processes 
will tend to freeze into a low-thermal-conductivity blan- 



ket that could allow Tq to exceed Ts . A blanket formed by 
freezing Earth's cu rrent atmosphere would on ly increase 
To by about 4 K (iLaughlin fc Adamsl [2000l ) . Isostatic 
adjustment (jFowleri Il990f r" however, will tend to cause 
some continents or islands to rise above the layers of wa- 
ter and ice on an Earth-like Steppenwolf planet with an 
active mantle, allowing volcanoes to continuously emit 
gasses that can freeze onto the water ice surface. Here, 
we will consider carbon dioxide because it is likely to 
be outgassed in significant quantities from an Earth-like 
planet and it supports a stable layer with relatively high 
base temperature relative to other common gases. To 
find an upper bound on the impact of a frozen gas layer, 
we will assume that the layer is Rayleigh-Taylor stable 
with respect to the underlying water-ice, so it remains 
as a surface blanket. 

To find the temperature at the carbon-dioxide layer 
base, we note that the thermal conductivity of car- 
bon dioxide again scales as T~^ (here, in a more lim- 
ited regime), but with constant of proportionality A « 
100 W m-i (jSumarokov et al.l[200l . We find that the 
maximum temperature supported is robustly w 220 K 
for Earth-like Steppenwolf planets and is determined by 
the weak temperature dependence of the melting curve 
(|Giordano et al.l 120061 ). Setting Tq = 220 K reduces the 
required steady-state water-ice thickness by an order of 
magnitude. 

At Earth mass, the temperature at the bottom of a 
layer of solid CO2 reaches the melting temperature of 
CO2 for a layer thickness of w 2 km, or « 3 x 10^ kg m~^. 
Venus' atmosphere has a partial pressure of CO2 of 
W90 bar (« 10^ kg m~^), which is roughly the vapor 
pressure of carbonate rocks at Venus' surface tempera- 
ture, implying that there may be more c arbon locked in 
rock in equilibrium with the atmosphere (jPierrehumbertl 
|2010[) . The store of carbon in carbonate rocks in Earth's 
interior is uncertain, but the continental crust is esti- 
mated to contain the equivalent of « 7 x 10^ kg m^^ 
CO2 and the mantle may contain 2-4 times this amount 
(|Zhang fc Z indler 1993). Therefore it appears reason- 
able to assume that a Steppenwolf planet could have a 
sufficient complement of CO2 to significantly elevate Tq. 

Since t he viscosity of ice dep ends strongly on tem- 
perature (jBarr fc Showmanll2009[ ). if ice convection were 
to occur on a Steppenwolf planet, it would occur only 
in the lower, warmer ice regions and would be c apped 
by a "stagnant" conducting lid (jSolomatovl |1995() . We 
calculate the ste ady-state ice thickness w hen convection 
occurs following iHussmann et al.l ()2006[ ). who assume a 
Newtonian rheology. We o utline the solutio n here, but 
the reader should consult IHussmann et al.l (|2006[ ) for 
more detail. We assume that convection occurs below 
temperature Tc (at higher temperature) and determine 
Tc by assuming that the viscosity is reduced by a factor 
7 = 10 over the convecting region, where the viscosity is 
given by the relation, r]{T) = 770 exp[l(Tm/T — 1)], where 
r/o = 10^^ Pa s, Tm is the melting temperature, and 
/ = 25. Assuming that Tm = Th, we can solve for Tc. A 
Nusselt-Rayleigh number scaling (Nu = aRa^) yields the 
thickness of the convecting region between temperature 
Tc and Th, 
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where we evaluate the viscosity at the mean tempera- 
ture of the convecting layer, T = ^{Tc + Th), k ~ 1-47 x 
10-6 m^ s-i, a = 1.56 x 10"^ K'^, p = 917 kg m-^, 
fc = 3.3 W m-i K-\ a = 0.12, and /3 = 0.3. Taking k 
to be constant in the convecting layer is a reasonable as- 
sumption given that the temperature is nearly constant 
within it. We add the thickness of the convective layer 
given by Equation ([2]) to the thickness of the stagnant lid, 

Hiid ~ y log (^if^: to find the total ice thickness when 

there is convection. When the total thickness with con- 
vection exceeds the conductive thickness given by Equa- 
tion ([T]), we use it for the ice thickness. This corresponds 
to a critical Rayleigh number of roughly 1000 within the 
convecting layer. 

Given that the ice composition on a Stcppenwolf planet 
is unknown and the ice material properties under appro- 
priate conditions are poorly-constrained, our convective 
calculation should be viewed as a rough estim ate. For 
examp le, following iHussmann et al.l ()2006( ) and iFu et al.l 
(|201C1|) , we have assumed Ne wtonian ice flow, which 
may or may not be realistic ()Barr fc ShowmanI [20091) . 
Creep mechanisms with stress-dependent viscosity, how- 
ever, should yield results rou ghly similar to Newtonian 
flow (jKirk fc Stevensonlll987t ). Another source of uncer- 
tainty is the appropriate ice grain size, the typical size of 
individual components of polycrystalline ice. In general, 
larger grain sizes correspond to higher viscosities, mak- 
ing convection less likely. The value of, 779 that we have 
used, 10^^ Pa s, corresponds to an ice grain size of wO.l 
mm (iBarr fc Showmai]||2009f) . For reference, the typical 

ice grain size on terrest rial i ce sheets is 0.1 10 mm 

(|Barr fc Showman|[2009f) and lFu et all (I2OIOI ) assume an 
ice grain size of 0.2 mm for their calculations of con- 
vection on icy extrasolar planets. Competing estimates 
of ice grain size on Europ a, however, place it at e ither 
0.02-0.06 mm or > 40 mm (|Barr fc Showman|[2009t) . Al- 
though the appropriate ice grain size on a Steppenwolf 
planet is a source of uncertainty, we have used a relatively 
small value, which is conservative in that it makes con- 
vection more likely to occur. Finally, the constant /S can 
take different values between 0.25 — 0. 33 depending on 
the g eometry and boundary conditions (jHussmann et af] 
l2006f ). We have used a fairly high value of /3, which leads 
to a conservatively strong scaling of i?conv with M . 

A Steppenwolf planet with mean ocean depth greater 
than the steady-state ice thickness, accounting for the 
ice-water density difference, will have an ocean under its 
ice layer. Figure [T] shows the required ice thickness to 
achieve liquid water in several scenarios. It is expected 
that the abundance of water on planetary surfaces varies 
greatly (jRaymond ct al. 2007,), but we can consider a 
simple scaling to understand how the conditions for liq- 
uid water scale with planet mass and water complement. 
At fixed planetary water mass fraction and fraction of 
water at the surface (rather than in the mantle), the 
depth of the ocean scales approximately as the typical 
depth at Earth mass (D*) times (M/Me)!-^^ Figure[T] 
shows contours for planets with water complements that 
are 0.1, 1, 10 and 100 times that of Earth. Combining 
this scaling with Equation ([T] and ^ , we find that if a 
Steppenwolf planet is similar to Earth in water mass frac- 
tion [D* = w 4 km), radionuclide composition, age, 
and has no frozen CO2 layer, it must be ~ 3.5 times more 
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Figure 1. Ice thickness that is required to support liquid water at 
its base as a function of planetary mass, assuming: (1) solid black: 
a melting point of 260 K and surface temperature required to ra- 
diate the geothermal heat flux in steady state, (2) dark solid gray: 
a melti ng point of 176 K at the euteetie mixture of water am- 
monia I IGrasset Hi SotinI 1199611 and surface temperature required 
to radiate the geothermal heat flux in steady-state, and (3) light 
solid gray: a melting point of 260 K and water-ice top tempera- 
ture along the melt curve of carbon dioxide at « 220 K (providing 
an upper bound on the possible impact of an insulating carbon 
dioxide layer) . The ice thickness is calculated using the conductive 
(Equation Q) or stagnant-lid convective (Equation |(2]|) solution, 
as appropriate. Convection increa ses the required thickne ss at low 
mass, and follows the method in IHussmann et al.l 1)20061 ). In the 
case where ammonia has suppressed the melting point by nearly 
100 K, we find a modest change in the conductive regime, but that 
convection is suppressed (at fixed water-ice rheology). The dashed 
black lines show the adjusted ocean depth (ocean depth multiplied 
by the ratio of the density of water to that of ice) as a function of 
planetary mass using the simple scaling M^~^'^ for water comple- 
ments that are 0.1 (lowest line), 1, 10, and 100 (highest line) times 
that of Earth, which is taken to be 4 km at 1 Mq. A subglacial 
liquid ocean is possible when the adjusted ocean depth exceeds the 
ice thickness. 



massive than Earth to sustain a subglacial liquid ocean. 
Contaminants which suppress the melting point have lit- 
tle effect on the conductive ice thickness, but significantly 
reduce the ability of the ice to convect. If a Stcppen- 
wolf planet has ten times more water {D* = lOD^) than 
Earth or if it has a thick frozen CO2 layer which reaches 
the maximum temperature of « 220 K at its base, the 
planet must be only w 0.3 times Earth's mass to have a 
liquid ocean. 

3. OBSERVATIONAL PROSPECTS 

We expect that detection of refiectcd sunlight in the op- 
tical wavebands and IR follow-up present the only viable 
observational choices in the near term . For a single-visit 
hmiting magnitude r « 24.7 of LSST (jJones et al.ll2009D 
and compara ble r = 24 in the nearer-term Pan-STARRS 
Jewittl[2003[) ). and albedo of 0.5, the limiting distance 
out to which an object can be detected with reflected 
sunlight is « 83O(r/i?0)^/^ AU. The Palomar survey of 
-12,000 deg^ to magnitude 21.3 (jSchwamb et all 120091) 
has discovered no such objects outside of the smaller 
trans-Neptunians such as Sedna. 

The baseline requirement to identify a Stcppenwolf 
planet is a detection of thermal emission in the far-IR. 
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The flux at the Wien maximum is 

where d is the Earth-object separation. At 10 Af®, 
Tg = 46 K so that Amax = HO ^m. Here, the Herschel 
PACS instrument reache s a 40 beams /sour ce confusion 
limit at a flux of «2 mJy (|Berta et al.l ()2010[) , suggesting 
a limiting dista nce of w4000 AU; PAC S reaches 10 mJy 
at 5cr in 1 hr (|Poglitsch et al.l l2010f )). Higher resolu- 
tion is required to progress to lower flux limits. At 
200 /im, the planned 25 m Cornell-Caltech Atacama Tele- 
scope (CCAT) would reach the source confusion limii|3 at 
0.36 mJy. 

Photometric microlensing has also been proposed as a 
method to detec t rogue planets throughout the galaxy 
(|Han et al.l 120041) . If each stellar system ejects one 
planet, a survey like the Galactic Exoplanet Survey Tele- 
scope could anticipate ^20 detections of rogue planets 
(|Bennett fc Rhiel[200llBennett et al.|[2010l ). Typical dis- 
tances to these objects would exceed the capabilities of 
follow-up that could elucidate their nature. Free-floatin g 
super- Jupiters have been discovered (jBihain et al.ll2009l ). 
but these represent a different class of objects. 

4. DISCUSSION AND CONCLUSION 

A Steppenwolf planet's lifetime will be limited by 
the decay of the geothermal heat flux, which is deter- 
mined by the half-life of its stock of radioisotopes (^°K, 
238^ 232rp]^^ and by the decay of its heat of formation. 

As these decay times are ^1 5 Gyr, its lifetime is 

comparable to planets i n the traditional ha bitable zone 
of main-sequence stars (jKasting et al.llT993l ). 

If a Steppenwolf planet harbors life, it could have origi- 
nated in a more benign era before ejection from the host 
star. Alternatively, after ejection, life could originate 
around hydrothermal vents, which a re a proposed loca- 
tion for the origin of hfe on Earth (|B aross fc HoffmanI 
Il985f l. If life can originate and survive on a Steppenwolf 
planet, it must be truly ubiquitous in the universe. 

We have shown that an Earth-like rogue planet drift- 
ing through interstellar space could harbor a subglacial 
liquid ocean despite its low emission temperature, and 
so might be considered habitable. Such an object could 
be detected and foUowed-up using current technology if 
it passed within 0(1000 AU) of Earth. 
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